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ABSTRACT 

We study the corrections to the conformal evolution of primordial magnetic 
fields after recombination, produced by the fall of the ionized fraction of matter 
into the dark matter gravitational wells. This effect enhances the field around 
the protostructures being formed, and might help to understand the fields 
observed in galaxy clusters and hydrogen clouds. 
subject headings: gravitation - magnetic fields - cosmology: theory. 
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1. Introduction 

The goal of this paper is to study the evolution of cosmic magnetic fields interacting 
with charged matter inmediately after recombination. At this time most ordinary matter 
combines in neutral hydrogen, but still a non negligible ionized fraction remains. On the 
other hand, as we shall show, the cosmic microwave background radiation (CMBR) is 
already too cold to provide sustantial dissipation, and so the ionized matter, as well as 
the neutral one, falls within the gravitational potential wells generated by dark matter 
overdensities. The infall of charged matter drags the field lines, thus resulting in a field 
amplification in and around the proto structures being formed. Thus a mechanism results 
which may provide a preamplification of a primordial magnetic field, before other best 
known processes, such as the galactic dynamo QZerdovich et al. 1990| ), become operational. 



The primeval motivation for the search of these corrections is to try to understand the 
presence of magnetic fields of intensity and coherence similar to the galactic ones, in galaxy 
clusters and much less evolved gaseous systems such as the damped Lyman-a clouds 
( Prcn fc Wolfe 1995Q . Although the existence and ultimate characteristics of those fields 
can hardly be explained by gravitational collapse alone, it is important to determine if 
this process may help to relax the requirements on the mechanisms for primordial field 
generation. 

The existence of relatively intense and coherent magnetic fields in galaxies, clusters of 
galaxies and hydrogen clouds is at present theme of intense research. There is a variety of 
proposals aimed at explaining their existence and characteristics, but none is completely 
successful. They can be divided in two main sets: one takes the point of view of local 
generation, and the other of the amplification of a primordially generated field ( [Zerdovicn 



et al. 19901 , |Harrison 1973; , |Pudritz fc Silk 1989| , |Howard fc Kulsrud 19961) . While both 



mechanisms compete on equal footing in explaining the presence of magnetic fields in low 
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redshift objects, like our own galaxies, the second approach seems more apropriate for 
highly redshifted, less developed systems, like the Lyman-a clouds. A model based on 
amplification mechanisms requires two stages, one to create the field and another to amplify 
it. For the first process there are several proposals (Polgov 1993| , Polgov fc Silk 19931 , [Ratraj 
19921 , [Mazzitelli fc Spedalieri 19951 , ICalzetta et al. 1998| , |Cheng fc Olinto 1994| , |Sigl et al 



TUWA [Lnqvist fc Ulesen 1994|, |Vachaspati TWT], [Hogan 1983;, |Martin fc Davies 1995|), some 



of which create fields strong enough but of coherence length much smaller than observed, 
while others generate coherent fields but of very weak strength. As for the second stage, 
the "turbulent dynamo" was a preferred amplifying mechanism for a long time flZerdovicE 
et al. 1990|) . But once again we must question it after the detection in damped Lyman-a 
clouds (|Oren fc Wolfe 1995] , [Wolfe et al. 1992|) of magnetic fields of the same intensity 
and coherence as the galactic ones. Another promising mechanism are inverse cascades of 
magnetic helicity ( Plesen 1997| Brandenburg et al. 1997[) , but, while this effect is quite 



well understood in ordinary magnetohydrodynamics ([Moffat 1961|) , its generalization to 
cosmology is still an open problem. 

It is well known from the theory of Large Scale Structure Formation, that galaxies are 
formed due to the infall of barionic matter into the gravitational potential wells of dark 
matter, the formation of those last beginning at the moment of equality between matter and 
radiation (z ~ 20000). Although they are not "deep" wells by the decoupling of barionic 
matter and radiation (z ~ 1100), as we know from the anisotropy of the CMBR, their 
existence cannot be disregarded. Despite this fact, current estimates in the literature often 
assume that, between decoupling and the starting of the dynamo mechanism or non linear 
gravitational collapse, their propagation takes place in a homogeneous, FRW Universe, 
and consequently that their evolution law is B oc Bq/o 2 , where a is the scale factor of the 
Universe and B Q an initial field intensity. In this work we study in more detail the evolution 
of a primordial magnetic field in a matter dominated, inhomogeneous Universe, with the 
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aim of estimating the corrections to this conformal evolution due to the evolution of the 
background perturbations. 

Long time ago Ya. B. Zel'dovich et al (1983) investigated the effects of a magnetic field 
over anisotropic self gravitational collapse of matter. They also estimated the amplification 
that a primordial field would undergo by that collapse, finding a factor of O (10 3 — 10 4 ). 
Working backwards from present estimates, the inferred intensity for a primordial field 
results then B ~ 1CT 9 — 1CT 10 Gauss, much larger than the ones predicted by the available 
mechanisms of primordial generation, for the required coherence scale. 

Our proposal differs from Zel'dovich's in that we do not study self gravitational 
collapse, but rather the infall of charged matter into preexistent potential wells created by 
dark matter. Consider a spatially flat Friedmann - Robertson - Walker Universe endowed 
with adiabatic perturbations. The main matter component is cold dark matter, which for 
our purposes only means that the smallest structures to be formed are galaxies. At the 
moment of equilibrium between matter and radiation, the matter dominated era of the 
Universe begins and the adiabatic perturbations in the dark matter component can grow, 
creating an inhomogeneous Universe. On the other side, the rest of the matter component, 
the baryons, remain coupled to photons until the temperature of the Universe is low enough 
to let electrons and protons combine. Then neutral hydrogen emerges while a small ionized 
fraction of O (1CT 5 — 1CT 7 ) remains. 

Any magnetic field present in the Universe will evolve coupled to the charged matter 
plasma. This means that before recombination, the field will be tied to the electron-photon 
plasma, and will remain coupled to the remaining ionized matter after that moment. Before 
recombination the plasma is hot enough to be considered relativistic and therefore its 
evolution can be followed without taking into consideration the background inhomogeneities. 
After recombination, however, the situation is the reverse: the temperature is too low to 
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neglect the effect of the perturbations and therefore we must study the evolution of the 
magnetic field coupled to the remaining ionized matter, in an inhomogeneous Universe. Of 
course the ionized plasma will be coupled to the CMBR as well as to any primordial field; 
one of our objectives is to show that, in the era under consideration, the interaction with 
the CMBR does not affect sustantially the collapse process, and that our plasma can be 
consider as a perfect fluid with infinite conductivity. 

Several authors have studied related physical situations. For example Ryu et al (1998) 
and Blasi et al (1999) studied the effect of inhomogeneities in the matter distribution of 
the Universe on the Faraday rotation of radio signals and light from quasars respectively. 
The inhomogeneities considered by the former authors are large scale filaments and sheets 
(generated by a hydrodynamical code), to which the magnetic field is glued. They assumed 
that a seed field is locally generated by a battery mechanism and further amplified by 
streaming and turbulent motions. Blasi et al, on the other side, considered redshifts for 
which the inhomogeneities are well described by the Ly-a forest. The magnetic field is also 
glued to this structure and they assume it scales with the electron density as B oc n^ 3 . In 
both works an upper limit for the magnetic field strength in the structures is derived. 

Our work differs substantially from the one by Blasi et al (1999) in that we do not 
assume a given relationship between magnetic field and ionized matter density, rather we 
derive it. The difference with the work by Ryu et al (1998) is that we do not assume a 
locally generated magnetic field, but a primordially generated one. Besides, we study the 
evolution of the magnetic field in a universe permeated by primordial inhomogeneities 
that are still in their linear regime of evolution, which means higher redshifts than in the 
mentioned works. 

We have therefore the following scenario: a weak primordial magnetic field is coupled 
to the matter that remains ionized after the recombination of hydrogen, and the whole 
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system propagates in an inhomogeneous universe. The matter component falls into the 
background overdensities due to gravitational attraction, dragging with it the magnetic 
field lines. This drag produces a local enhancement of the magnetic field in the overdense 
regions, deviating from the evolution to be expected in a pure FRW Universe. Our purpose 
is to calculate this enhancement and to see if it is intense enough to explain (up to certain 
level) the detected field intensities. 

The paper is organized as follows: In section 2 we describe the geometry of the 
background Universe, in section 3 we describe the matter component the magnetic field is 
tied to, discussing in some detail the possible interactions with the CMBR. In section 4 we 
study the evolution of the magnetic field and summarize our main conclusions in Section 5. 



2.1. The geometry of the background 

Consider a matter dominated flat FRW Universe, endowed with adiabatic perturbations. 
In the longitudinal gauge, the metric tensor for this Universe reads (Prandenberger et al. 



where a (rj) = (rj/r]o) 2 is the scale factor, rj the conformal time and 7^ the spatial flat three 
metric, 7^ = diag(l, 1, 1). We have normalized the conformal factor to unity at the time of 
decoupling, given by 77 = r] ~ OAh^ 1 x lO^GeV^ 1 . Greek indices run from to 3, while 
latin ones take the values 1, 2, 3. The function ip depends on r and rj and can be interpreted 
as a generalized newtonian potential. Its amplitude can be determined by the level of 
anisotropies in the background radiation caused by the Sachs- Wolfe effect at decoupling. If 



2. Baryon free fall 



1992D 




(1) 
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ib describes a growing perturbation then it must depend only on r ( [Brandcnbcrgcr et al 



1992| ). In the linear regime, this mode will grow as ~ a(r])ip(r). 



From this expression we can find the four velocity of the background fluid flow and 
construct a 3 + 1 description. This is needed because we will consider the presence of a 
magnetic field, and its definition relays on the specification of the spatial surfaces of the 
spacetime. 

In a coordinate system where peculiar velocities are zero we have 

TTII dx» (I -ib J\ 

This is normalized by 

= -1 (3) 
Finally we build a projection tensor for the spatial three surfaces as 

hT = + UW (4) 

As this fluid has zero vorticity, the spatial three surfaces are orthogonal to the fluid 
four velocity QEjEi fc van Elst 1998| ) 



2.2. The plasma component of the Universe 

We want to study the evolution in the matter dominated era of the Universe of a 
primordially generated magnetic field, considering that in that period of time the process 
of structure formation is taking place as well. At the beginning of matter dominance 
(z ~ 20000), the background perturbations begin to grow, while the baryons remain 
coupled to the cosmic background radiation. When the temperature of the Universe is 
low enough, electrons and protons recombine (z ~ 1100) and most of the barionic matter 
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becomes electrically neutral. It can then start falling into the already formed dark matter 
overdensities, which are described (in the linear regime) by the function if>(r). There 
remains a ionized fraction of barionic matter of order ~ 0(1O~ 5 — 10~ 7 ). 

We assume that any magnetic field present in the Universe interacts only with the 
barionic component of the matter. With the electrically neutral fraction, the interaction is 
very weak and is described by B.m where m is the atomic magnetic moment. With the 
ionized fraction, the interaction is given by the Lorentz force, being therefore stronger than 
the former. We will study the system of ionized barionic matter coupled to a preexisting 
magnetic field, in the matter dominated epoch, by means of the equation 

T£ = (5) 

where is the stress energy tensor for the matter and T^ 1 is the Lorentz force density, 
exerted by the magnetic field on the charged particles. 

We shall investigate the effects on the preexisting magnetic field of the free fall of 
matter into the newtonian potential wells described by if) (x) . To this end, we shall consider 
a geometry where if) = if> (z), and correspondingly the motion of matter is along the z 
direction only. In this first Section, moreover, we shall make two simplifying assumptions, 
namely, we shall consider the plasma as a perfect fluid, and we shall neglect the right hand 
side of Eq. fl5|). We shall show the validity of the first approximation towards the end of 
this Section; the second one will put limits on how far the model can be relied upon as a 
description of collapse, and will be discussed later on. 

We shall treat the energy density Se, pressure Sp and four velocity Su^ of the ionized 
plasma as perturbations on the corresponding quantities for the dark matter, namely 
e = e + Se, P = p + Sp, and W 1 = + Su^ 1 . We shall call SAf any conserved quantity proper 
of the ionized fraction, for example the leptonic or barionic numbers, and we will consider 
that this quantity does not vary with time, i.e. the ionisation fraction remains constant 
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( |Lucchin fc Coles 1995 , Zerdovich fc Novikov 1983 ) Since the gravitational effect of this 



perturbation is negligible, we shall not consider the perturbations to the geometry that it 
could produce. We are therefore describing the free fall of the plasma into the potential 
wells created by the dark matter. The full energy - momentum tensor for the plasma is 

= + 8Tf (6) 

where 

= eUW + -pbT (7) 

STj? = 5eU^U v + {e + p) {5u^U" + U^5u u ) + Sph^ (8) 

As stated before, we have 

xr = o (9) 

We will also need the equation for the conservation of the number of matter particles, 
namely 

Z£ = o (io) 

where 

= SAT (U^ + Su") (11) 



2.3. Timelike and spacelike projections of Tf^ 

Our ultimate goal is to study the evolution of the matter component and the magnetic 
field under the gravitational perturbations created by dark matter. As the specification 
of what magnetic and electric fields are is tied to the definition of spacelike surfaces and 
their corresponding timelike curves, we must work in a 3 + 1 formalism. For this purpose 
we project equation ([]) onto the spatial surfaces, by h% (eq.(f|)) and along the worldlines 
given by f/ M (eq. (0)). Performing the covariant derivatives with the Christoffel symbols 
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corresponding to the metric tensor Eq. (Jl]) (see Appendix) and replacing the explicit 
expressions for the four velocity we obtain the conservation laws to leading and first order 
as follows: 

zeroth order 

e + 3-(e+p) = (12) 
a 

For the dark matter component, moreover, p = 0, and we get the expected behavior 
e ~ cr 3 . 

first order 

(e + p)- 3^ (Se + Sp) -(e+p) 6u J }j = (13) 

tl( £ + p ) + ^ + (e + p) — + (e + p)(-5u i + 5^-5u i ) =0 (14) 
a z a z a la a z J 

We now replace the zeroth order equation and write 5e = pe, 5u^ = /a and p = c^e, 
where c s is the speed of the sound in the medium. We will also neglect the terms with tp as 
they correspond to decaying modes of the background perturbations, leading to 

p+ (1 + ^)^=0 (15) 

V + j^rf + - 3c ') -/ + ^ = (16) 

In what follows, we neglect the effect of the pressure. This is justified in view of the 
low temperature of the matter. 



ib 1 r 

— e be + 3 

a a 



a z a 
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2.4. 



Particle number conservation 



Developing the covariant derivative of eq. (|T(]) we obtain 



8Af (1 - V>) + 35JV- (1 - V) + (M/V) ;J -2AA/V 




(17) 



Taking <5jV = 5n/a 3 the previous equation reads 



(1 — iji) 8h + u 3 5n ,j —25nilj,j u 3 = 



(18) 



2.5. Study of the evolution of p and u J for a dustlike gas 



We are now ready to analyze the effects that the background perturbations may have 
on the density and velocity of the plasma. Current numerical investigations of large scale 
structure formation show that gravitational colapse of matter is highly anisotropic, giving 
rise to a complicated network of rich and poor clusters connected by filaments with void 
regions in between ( Pond et al 199(f ). This is a difficult process to study analytically so 
we will retain one simple feature, namely that gravitational collapse is not symmetric, but 
takes place mainly along a preferred direction, the resulting structures being pancake or 
cigar-like shaped. This fact allows us to introduce another simplifying hypothesis: we will 
consider that all quantities depend only on the coordinate along which gravity acts, let it be 
the z coordinate. This simplification contains a very important physical assumption about 
the perturbations: they are vorticity free. Were we to consider two or three dimensional 
gravitational collapse, vorticity must be taken into account. 

When we make the mentioned approximations we are left with 







(19) 



+ -u i + u i = 



(20) 



a 



13 



We write u = (u x , u y , w 2 )and obtain 



p + d z u z = 



a a 
u x H — u x = u v H — u v = 
a a 

CW + M 2 + = 

a 



(21) 

(22) 
(23) 



For the initial conditions, u(z,r] ) = 0, and p(z,r) ) = p Q , the solutions to equations 

(24) 



T|) and (H) read 



p = po + fkvw 



'-V+ 


% _ 




77 






[(?)■ 


_ ^ 






= 





(25) 
(26) 



These are the equations which describe the free fall of the plasma into the potential wells. 



2.6. Discussion 

Equations fll8|) and (p0| ) are the goal of this Section, namely, they describe the evolution 
of the number density and velocity of the plasma as it falls on the dark matter overdensities. 
Since the magnetic field lines are tied to the charges, we expect these will be dragged along 
with matter, producing an amplification of the field in the overdense regions. The study of 
this effect is the subject of the rest of this paper. 

However, before we go on it is right that we discuss the first of the approximations 
under which eqs. (|i~9l) and ( ^0|) were derived, namely, the assumption that the plasma 
could be described as a perfect fluid. In reality, the plasma interacts not only with the 
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gravitational field of the dark matter, but also with the photons in the CMBR, and this 
interaction could in principle result in dissipative behavior. 

In the epoch we are interested in, the plasma and the CMBR are not in equilibrium, and 
therefore the usual estimates of the dissipative effects induced on the plasma do not apply 
( |Weinberg 1971| ). Indeed, the mean free path for the photons is l m f p = (n c cr T ) _1 , where n c 
is the number density of the electrons in the ionized fraction of matter and gt = 87re 2 /3m^ 
is the Thompson cross section, (due to the low temperature of the background photons, the 
main interactive process is Thompson scattering). Their values are oj> — 1.2 x 10 5 GeV~ 2 , 
n c ~ 0.4 x 10~ 45 a(^)~ 3 Ge^ 3 and therefore l m f p ~ 2 x lO^a^^GeV" 1 . On the other side, 
the size of the particle horizon is dh = rf /rfc. We estimate the number of collisions that a 
photon suffers per Hubble time as n co u = dh/l m f p ~ 20 a(ri)~ 3 / 2 , which is too low to sustain 
equilibrium. 

Still, the plasma must feel some kind of friction as it moves accross the CMBR. 
Since the plasma is electrically neutral, any baryon motion is matched by a corresponding 
electron flux. The Thompson cross section is 10~ 6 orders of magnitude smaller for protons 
than for electrons, then the interactions of photons with electrons dominate, and we can 
consider that the electrons suffer a friction due to the photons. Therefore, instead of the 
conservation equation @ we should write 



T£ = C; 2£ = (27) 

where X M is the particle four current. We also have the corresponding equation for the 
photons 

T% = -C (28) 
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From the first principle of thermodynamics we have 

S" = - (3 CU T^ U - (3 1V T^ - clP (29) 

where is the entropy flux, (3g = U^/T c , (5^ = U^/T^ , the subindex c corresponding 
to the "charged" matter and 7 to the CMBR. $ M is the thermodynamic potential (Israel 
1988), whose derivatives are 



<9$M QfoH <9<T)M 

d(3 cu ' <9/3 7i/ 7 ' 9a 



By the second law of thermodynamics we have 

S£ = (/V-«C>o (3i) 

which is satisfied only if 



= Ctf* (-/% - + C 2 - ft) 



(32) 



7 

where /?? = — /3 7At /3ijf and C\,Gi > 0. In terms of the four velocities and temperatures eq. 
(j32|) reads 

C = -^TTjT [T e + TJJJT] -^[U» + UWU V ] (33) 
We see that the first term in the r.h.s. of eq. (|33| ) can be interpreted as the heat interchange 



between the photons and the plasma, while the second one as the momentum transfer 
between them. 

To estimate the C\ and C2 coefficients, let us place ourselves in the rest frame of the 
CMBR. Since in this frame the photon bath is isotropic, the bombardment of electrons "at 
rest" by photons averages out, and any net effect is solely due to the motion of the electrons. 
Each time an electron strikes a photon, the later gains a momentum Ap 7 ~ 1/A ~ T 7 . The 
electron loses a momentum Ap e ~ — T 7 and changes its energy by AU = — t> e T 7 , where v e is 
the velocity of the electrons. The number of collissions per unit time that an electron suffers 
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is dn = n^aTVg. If the electron number density is given by n e , then the net force density 
excerted by the photons on the electrons is given by F^ oll ~ n e n 7 o"rT 7 f e , and the mean 
energy loss is Q = rien^aiT^i^. Performing the average over all electrons, and comparing 
with eq. (|3~3"D , we conclude 



C\ n c n 7 a T T 7 



(34) 



C 2 

— = n c n^a T T^ (35) 

J- c 

We should now project eq. ( |3~3"D along the fluid flow lines and onto the orthogonal 
spatial surfaces and replace the expressions in the r.h.s. of eqs. (|X3|) and (|14"D respectively. 



For the timelike projection only the first term of eq. ([33|) contributes while for the spacelike 
projection we need the second term of that equation. Considering e as the critical density 
of the Universe, using also that e cr u = e^\ t /a 3 (r]), n = n (°Va 3 (^) and that T = T^/a(rj), 
where with the supraindex (0) we refer to quantities at the epoch of recombination, we have 
that the equations ( |T§| ) and ( PU|) now read 



™ e 4rU 3 (f?) 



P + = - ° Jo) l,\ ^ ~ T i) (36) 



t + -u l + u l = c —^- ^— u l 37 



Replacing the figures, ~ 10" 38 GeV 4 nf ~ 1(T 33 GeV 3 , a T ~ 10 5 GeV~ 2 , T<® ~ 10 
GeV, we obtain for the r.h.s. of eq. (|36| ) 
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1,(0)^(0)^^(0)2 in _49 

7 \ T 7 * -TT^GeV (38) 

where we have assumed T c <C T 7 . We compare this quantity with the divergence of 
the velocity in the l.h.s. of eq. (|T9|), using eq. (p5| ). Asuming <9 2 ~ Zq 1 where z is a 
characteristic scale of the background inhomogeneity, which for a galaxy is z^ ~ 10 35 
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GeV -1 , we have T]q/Zq ~ 10 2 . Recalling that initially the peculiar velocity is zero, we have 
that the effect of the heat transfer becomes negligible inmediately after the matter starts 
falling into the dark matter overdensities. 

For equation (|20|) we compare the second term in the l.h.s. with the r.h.s.. For the first 
we have that the coefficient of u l reads a(rj)/a{rj) = 2/rj = 2r]Q 1 a~ 1 ^ 2 (r]). For the coefficient 
of u l in the r.h.s. we have 



,(°) 

'crif- 



This factor will be smaller than the corresponding one in the l.h.s. of eq. ( E^) when 
a{rff' 2 > 1, and this constraint is always satisfied, because we have a(rj) > 1. We therefore 
have that in the whole time interval considered in the paper, the equations for a perfect 
fluid apply. 



3. Evolution equation for the magnetic field 

As we have seen in the previous Section, after decoupling the partially ionized plasma 
falls into the newtonian potential wells created by dark matter. In its infall, it will dragg 
the magnetic field lines, producing an enhancement of the field intensity, proportional to 
the depth of the potential well. Our goal in this Section is to study this effect. In the last 
subsection, we will discuss the limitations on our analysis. 



3.1. Maxwell Equations 

Let us begin by deriving the Maxwell Equations that satisfy the electromagnetic field 
in the spacetime we are working in. These were put in the 3 + 1 form for the first time by 
G. F. R. Ellis (1973) and latter by K. Thorne & D. MacDonald (1982). 
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We begin by denning the field strength tensor as usual 

Ffj, v = d„A u - d v A„ (40) 

where A m is the electromagnetic four potential. The electric and magnetic fields are given 
by 

= F^U V = -F^U U (41) 
= fi^lTF** (42) 

We can rewrite the field strength tensor in terms of these fields as 

F, u = U,£„ - S,U U - Tj^pirBP (43) 

and with the indices up 

where is the four velocity of fiducial observers, which in our case is given by equation 
(@). The Maxwell equations can be written in covariant form as 

F^ v = J» (45) 
V^F^ = (46) 

We will need the projections of these equations onto spatial surfaces and along the 
four velocity. We also transform the electromagnetic field and current as £ l — > E l /a 2 , 
B l — > B l /a 2 , J l — > J l /a 3 . The previous set of equations then reads 

-E k , k +3^, k E k = -(l+^)J° (47) 
- E* (1 - + rf j m \{l - if,) B\ 3 -to/,,; B k ] = T (48) 
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r>' ; ■■ r> J c, 



(49) 



V\ 



1 _ ^) E k ,j -3^ £ fc l + (1 - V>) ^ = 



These are the equations that we will use in the main part of this work. 



3.2. Plasma conductivity and the magnetic field 



(50) 



The currents and fields appearing in Maxwells equations above are further related by 
Ohm's law 

J" + WWJ V = aF^Uu (51) 

where is the fluid four velocity given by = U^ 1 + 5u^. Keeping only first order terms 
we have 



T - u l p e = a E l + if\ r vPB r 



pr 



(52) 



In our case, it is appropriate to consider the plasma as perfectly conducting. This 
is essentially due to the fact that any charge separation would lead to an electrostatic 
attraction much stronger than the gravitational forces considered so far, and thus it cannot 
be sustained over any macroscopic lapse. For an infinite conductivity, the electric field 
in the rest frame of the plasma must vanish, and this yields an extra relation among the 
electric and magnetic fields, and the plasma velocity in that frame, namely 



E i = -r] i0 ^u p B r 



pr 



(53) 



For equation (50) we have 



;i - ip) V x (u x B) - 3V^ x (u x B) = (1 - tp) 



OB 

dr] 



(54) 
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This is the evolution equation for the magnetic field. 

Besides the arguments above, it is possible to give a direct estimate of the plasma 
conductivity that confirms it may be considered as infinite. Suppose an electric field were 
imposed on the plasma, generating an electron flux with velocity v e and a current j = en e v e . 
The field yields a Joule power P = Ej, which must be equal to the power dissipated into 
the CMB by the friction force (recall the discussion in the previous Section) 

P = 2n e n^a T v 2 e T = 7 x9 e (55) 

a(r}) 



Therefore 



from where we read the conductivity 



Ve = (56) 

Tricar 



o c = (57) 

Estimating its value with the same figures as before we obtain 

a c ~ 10- n a (rj) GeV (58) 

In order to determine if this conductivity is large enough to be considered as infinite, 
we may compare the conduction current with the displacement current generated by the 
collapse, d v E ~ E/i] , where t/q 1 ~ 10~ 38 GeV. We see that we can consider the conductivity 
of the plasma as infinite for all practical purposes. 



4. Evaluation of the magnetic field 

The evolution equation for the magnetic field can be simplified as follows: Write 
u = (0, 0, u z ) and B = B (77, z) . Then 

VxB = UB S) 9A,0) (59) 
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u x B = (-u z B y , u z B x , 0) 

V x (u x B) = {-d z [u z B x ] , -d z [u z B y ] , 0) 

Vip x (u x B) = u z B x , -ip, z u z B y , 0) 

We see that B z is not affected by the collapse. For B x and B y we can rewrite equation ([54] 
as (ommiting the subindex) 



(60) 
(61) 
(62) 



l-i)fB +d z (l-tpfBu 







(63) 



which expresses the conservation of the quantity (1 — ip) 3 B. 



At this point it is important to note that we cannot recast equation (|18[), which 
expresses the conservation of the number density that characterizes the ionized fraction, in 
a form similar to eq. fl6"3"|). This means that the relation between Sn and B might not be a 
simple one. Nevertheless since ip is a small deviation from inhomogeneity, we can obtain a 
conservation law from eq. flIB|) if we neglect the last term and the newtonian potential in 
the first term. We obtain 

Sh + (u z 5n) , z = (64) 

where 



u 



(65) 



Vo \ _ V_ 
K V ) Vo 

with f(z) = r] d z ip(z)/3. Performing the same approximations in equation (|63]), this 
equation reads 

B + [Bu z ] , z =0 (66) 

Since both B and Sn are assumed to be homogeneous at rjo, from eqs. ( |B"4[ ) and Q5"6] ) we 
may conclude that B oc <5n, at all times, where the constant of proportionality is given by 
the value of B/5n at n = r] , i.e. we have 

B 



B(rj,z) 



5n t 



-Sn (n, z) 



(67) 
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Observe that this is not the behaviour expected for a homogeneous universe, where w J = 
and consequently the solution to eq. (|66|) would be B = const. The solution to eq. (|B~4]) 
reads 



5n = -j^tH 

f{z) 



dz' ( 1 rf rjo 



f(z') +r]0 [2r 1 i + r, 2 / 
where H (C)is an arbitrary function obtained from the value of 5n at rj = rjo- The hypothesis 
that the magnetic field Bq is initially uniform over the horizon size, although simplistic, 
is valid for weak magnetic fields ( jCalzetta et al. 1998| ) and also helps to simplify the 



(68) 



calculations. The real initial large scale structure of magnetic fields is still under study 
Jedamzik et al (1998) have shown that linear magnetohydrodynamic modes suffer damping 
during recombination and neutrino decoupling, while Subramaninan and Barrow (1998) 
analyzed the non-linear case. They both find scales below which those modes are dissipated. 
Nevertheless their treatment refers to perturbations to a background magnetic field, while 
here we aim to study the evolution of that background field in an inhomogeneous universe. 

In order to achieve our purpose, we only need the newtonian potential profile ip(z) and 
to our ends, it is enough to consider that the dark matter is still in the linear regime (see 
below), the shape of the potential being then essentially given by the theory of primordial 
density generation ( Pardeen et al. 1983| , Guth fc Pi 1982 , Starobinskii 198% , [Hawking 



|1982|) . However, this theory does not yield a deterministic prediction, but only the relative 
probabilities of different density contrast profiles. Therefore we shall restrict ourselves to 
considering a few simple situations, namely, a power law density profile, a uniform slab of 
finite height, and a harmonic density contrast. 



Power law density contrast Let us consider dark matter density profiles of the form 



, £ [— Y 2 for Z < Zn 
8{z) = { U ' (69) 

for z > z 
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(70) 



They give rise to power law (generalized) Newtonian potentials, i.e.: 

a (£f f° r z < z ° 
a {pf o -p + l) for z > z 

where a has dimensions of length" 1 and its amplitude is a ~ 10~ 5 ,.and z determines 

the spatial extension of the perturbation, which can be the one of a galaxy 

(z ~ 1.42 x lO^GeV" 1 ), or cluster of galaxies (z ~ lO^GeV" 1 ) for example. 

Physically meaningful profiles require 1 < p < 2; p < 1 implies an infinite force towards the 

center of the structure and p > 2 means that the density increases towards the outer edges 

of it. The function f(z) reads 

for Z < Z 
for Z > Z 



m = 



PVpa 
3z 



(71) 



and 



dz' 
W) 



Z Q 



pr] a (2-p) 



i-vz 2 p for Z < Z 



34 

pr/ooc 



Z 



for z > z 



(72) 



As stated above, we evaluate the function H(Q by fixing the functional form of Sn (z, w) 
at i] = i] . In view of the lack of any prescription for it, we also choose it as a constant, i.e. 
Sn (z, f] ) = 5no We then have 



H[Q=5n 



34 



(2 - p)pr] a 
3z p n 



c 



(p-l)/(2-p) 



3zl 



pr] a (2 - p) 



(73) 



H[Q = 5n 



prioa 
32 ' 



otherwise 



(74) 



Finally the density profile 5n reads 



For z > z 



H[z,n] = 5n ^ (75) 



-24- 



f[z 



5n = 5n 



(76) 
(77) 



For z 2 - p - (a + £ 



3 



H[z,r}\ = 6n — — ; 

P-1 



o 



(78) 
(79) 

(80) 



For0<^<, 2 --^|^(^ + ^-|) 

prj a 



H[z, V ] 

M 



5m 



PVoot z 



3z 

■p-i 



2 (2-p)pa 2 ho IV 



34 



2 ^ 



<5n = 5n 



3 
1 



z: 



ZP~ 



o 



2 _ p (2-p)pa 2 (vo Iff 



(p-l)/(2-p) 



3^ 



3 



(p-l)/(2-p) 



51) 

«) 



We can see from eq. ( B7f ) that the magnetic field will grow with time for all exponents in 
the physical range 1 < p < 2. It remains unperturbed outside the dark matter distribution, 
and inside it freezes for long times with a profile given by eq. (0). For p = 1 we obtain 
B = Bq for the whole range of z. For example, for p = 3/2 and the density profile grows like 



5n 



1 f aril 

dn ^= < y/z H m 

'z I 4^ 7 



% ir _ 3 

J7 27/2 2 



(84) 



/or < y/z < sJzq~ 



aril 



4z 



3/2 



^0 



l77_ 2 

2% 2 



5n = 5n \ — 
V 2 



/or v^o ~ 



42 ; 



3/2 



Vo }_V_ _ 3 
V 2 % 2 2, 



< v 7 ^ < v^o 



55) 
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and 



5n = <5n ; for z > z 



(86) 



We now analyze the case p = 2. 



Homogeneous density contrast. Let us consider a dark matter density contrast profile 
given by 

£ for z < z 
for z > z 



5{z) = 



The generalized Newtonian potential can be written as 

a(£f for z<z 
2(2^-1) for z>z 



and 



The function f(z) then reads 

m 

dz' 



2ripa z 
3 ^ 

3zo 



/or 2; < 2 
for z> z 



3z. 



In (z) for z < z 

: /or 2 > 2 



2r? o 
3zp 



With the same initial conditions as before, the function H (() is given by 

if C< 



r 2ar] B f2r) a 1 

# [C] = <fa q 2 ex P i ^i2-C ; 



3*g 
iJ[C] = 5n 



{ 3^ 
2m]o_ 

3^0 



3z 2 

J_°_ln r 

o m ^o 

2r/ a 



otherwise 



The density contrast reads 



<5n = 5no exp 



2t7o" Ho . 1 ^ 2 



3z 2 \r, + 24 2 / 
2 



for < z < z exp 



2 ^ 

3z 2 V^? ' 2r? 2 2 / 



Ho 1 V 



(87) 



(88) 



(89) 



(90) 



(91) 
(92) 



(93) 
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Sn = 



r z 
on — ; 

z 

for z exp 



(94) 



3z 2 U 2r/ 2 



and 



<5n = 5n ; for z > z 



(95) 



Harmonic density profile Assume that we have a density distribution whose newtonian 
potential is given by 



1>(z) 









a 




K)] 









(96) 



This function is defined over the whole particle horizon, z is a characteristic scale 
of the background perturbations, that could be interpreted as the distance between two 
adjacent galaxies. The function f(z) reads 



f2irari \ . (Itxz 



and 



3^o r 
iTrar/o J 



dz 



3*3 



„— — lntan(7r — 

27iai] J sm{27iz/z ) 4ir z ari \ z 



(97) 



(98) 



According to what we did before, we have to fix the function H [(} by specifying the 
density contrast at r] . We therefore have 



Sn = 



3z 



27rm7o sin (^f) 



H 



3zi . ( z 
— In tan 7r — 

An^arjo \ z 



from where we have 



H [C] = 5n 



iTrarjo 


exp 




3z 


1 + exp 


z H 



(99) 



(100) 
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And the density profile is 

5n = Sn 



1 + tan 2 








)} 


1 + tan 2 




«p{^(¥ + iS-i) 


} 



(101) 



We find once again that the density contrast freezes for long times. This time, the final 
shape is given by 



Sn 



5n 



(102) 



sin2 

Obviously, the smaller the value of the sin function, the longer the time it takes to reach 
the value given by equation ( |102| ) 



4.1. Interval of validity of the calculations 

In the case of the power law density profile we see that the magnetic field diverges for 
z — > for all times. This is due to the fact that the density profiles diverges at that point. 
We will therefore analyze the limits of validity of the calculations for the case of uniform 
dark matter density contrast and of harmonic density. 

The first point is to specify until which moment the background perturbations can be 
considered in their linear regime. With the value a ~ 10~ 5 fixed by COBE, we have that 
the galactic scale entered its nonlinear regime at a redshift ~ 5, which corresponds to 
i] ~ 15?7o, where 770 — 0.4/i _1 x lO^GeV" 1 is the conformal time at decoupling. For a galaxy 
cluster, the dark matter is still in its linear regime (Peebles 1993). 

The other assumption we are going to check is neglecting the effect of the magnetic 
field on the evolution of the perturbations, namely the the r.h.s in equation @. The field 
affects the motion of the plasma through the Lorentz force 

f m = J n F mn (103) 
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F mn is the e.m. field strength tensor and J n the induced electric four current. Since F mn 
is spacelike to lowest order, the projection along the four velocity is negligible, and the 
orthogonal proyection, written in terms of the conformal magnetic field, yields [J x B]\ 
The current is given by Maxwell's equations fl47|) and fl48|) ; after substituting the electric 
field eq. (|53|), we obtain J° = (thus no charge separation) and 

J = d n u x B + u x d v B + (1 - tjj) V x B - 3V^ x B (104) 

(it is interesting to note that J is actually orthogonal to the macroscopic mass flow). For 
the symmetry of our problem, and neglecting pressure effects, the equation of motion for 
matter is changed into 



Byd z B y + B x d z B x 



d z ifj + ii z + -u z = 
a 

We have to find at what value of r\ the inequality 

Byd z B y + B x d z B x 



(105) 



\dM » 



CLEq 



(106) 



breaks down. Let us consider the uniform slab. Using equations (|8q ) and ( p4|) we obtain 



Zn 



> 



B 2 z 4 



6 z"e Q \ 7] ) \7] / 2az A e 

The correct bound is found by replacing the smallest value of the z/z coordinate from 
equation (|93|) to obtain 



(107) 



4 z 3 e \7}J \rj 0j 



> 



2aeo 



exp 



iaVo (Vo , 1 V 2 



3z 2 V 77 + 2 V 2 2 / 



(108) 



As stated before, is the critical density of the Universe at recombination, 
So ~ 8h 2 x 10 -38 GeV 4 . We take the value of the primordial field at decoupling as being 
given by the mechanisms proposed in the literature for a coherent field at galactic scale 
( Polgov 1993| , [Ratra 1992| , [Mazzitelli fc Spedalieri 1995| , |Calzetta et al. 1998| , |Uheng & 



-29- 



Olinto 1994, |Sigl et al. l"99?j |Knqvist fc Ulesen 1994!, |Vachaspati 199Tj polgov & 53E 



19931 , IHogan 1983| , |Martin fc Davies 19951 ), i.e. B today ~ 10- M Gauss ~ 10- 45 GeV 2 . We 
estimate its value at decoupling as B ~ (1 + Zdec) 2 B to day ~ lO^Gel/ 2 and we therefore 
have B 2 /eo<y ~ 10~ 36 . Replacing in equation (|108| ) we can check that the inequality is 



satisfied for i]/i] < 2.22 (Z ~ 200) for a galactic scale ~ 10 35 Gev _1 obtaining at that 
moment an amplification factor of the order 10 4 . For a cluster scale Zq ~ 10 36 Get> _1 
and equation ( |9"3"D is valid up to 77/770 < 17.5 for the same value of a, obtaining a similar 
amplification factor. 

It is also inmediate to check that the neglecting of the last term in eq. ( |TSD is a valid 
approximation for all times. 

We therefore conclude that the validity of neglecting the backreation of the magnetic 
field on the evolution of the perturbation fixes the time interval of validity of our 
calculations. 



5. Final Remarks 

We have studied the evolution of a preexisting cosmological magnetic field after 
recombination of hydrogen, in a background Universe described by a FRW geometry plus 
adiabatic perturbations, and considered that the field is coupled to the remaining ionized 
fraction. This matter falls into the background overdensities, dragging with it the magnetic 
field. This dragging results in a deviation of the behaviour of the magnetic field from 
the law B oc 5Af 2 ^ 3 where 5Af is a conserved number density of the plasma, even when 
the background perturbations are in their linear regime of evolution, as can be seen from 
equation ( |67D together with (p4|), (|93|) and ( |101| ). The background density profiles used, 
although simplistic, retain an important physical characteristic, namely the asymmetry of 
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the gravitational collapse. In the case of the harmonic potential, we could also think that 
it describes an alternating pattern of walls and voids ( Proadhurst et al. 1990| ), despite 
the difference with the scales used in this work. The rate of growth of the field depends 
on the characteristics of the background inhomogeneities and on the initial distribution of 
magnetic field intensities. The hypothesis of homogeneity of the field over the horizon scale 
at recombination can be sustained provided that the intensity of the field is sufficiently 
low, a fact achieved in most of the mechanisms of primordial field generation ( |Dolgov| 
19931 IKatra 1992| , |Mazzitelli fc Spedalieri 19951 , jUalzetta et al. 1998) , |Cheng fc Ulmto 



19941 |Sigl et al. 1997) , |Enqvist fc Olesen 1994! , |Vachaspati 1991] , Polgov fc Silk 1993| , |Hogan 



1983| , pVlartin fc Davies 1995|) : the larger the coherent scale, the weaker the field intensity. 



The level of amplification depends mainly on the scale of the background inhomogeneity, as 
stated in the last subsection, and have its maximum at the center of the structure. For the 
x — y plane, the initial symmetry remains because gravitational collapse is along the z axis 
only. As stated in the previous subsection, our calculations are valid up to rj ~ 2.22^0 for 
galaxies and for the whole time interval for cluster scales. 

Our main result is that it is not correct to consider the magnetic field evolves as 
B oc Bo/a(r]) 2 after recombination. We found that even in the linear regime of gravitational 
collapse, there can be nontrivial corrections to that behaviour. Even if the growth law were 
a power law, we could have that the field remains constant between recombination and the 
beginning of an amplifying mechanism, over the scale of interest. Clearly we do not claim 
to have solved the problem of the origin of magnetic fields in high redshifted systems, like 
the damped Lymann-a clouds, but if the scenario of primordial field generation plus further 
amplification is accepted, our results might help to relax the constraints on mechanisms 
for primordial field generation, because now the input fields for any further amplifying 
mechanism could be orders of magnitude stronger than the values that those fields would 
have if they had simply decayed as a~ 2 (r]). For instance, it is stated ( tZel'dovich et al. 1990| ) 
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that the present intensity of primordial magnetic field that can seed a galactic dynamo is 
B ~ 10~ 21 Gauss over a scale of 1 Mpc. Most of the proposed mechanism for creation of the 
field give a present value B ~ O (10~ 24 ), this means that with the usual conformal evolution 
of B, the value of the field at recombination would be B rec ~ 1CT 18 Gauss, and at rj ~ 2r/o, 
B ~ 10~ 19 , while with our calculations, the intensity of the field at that moment would be 
B ~ 1(T 16 Gauss. 

The analysis in this paper may be improved in several particulars, the most important 
ones being the inclusion of vorticity in the plasma and the extension of our calculations into 
the nonlinear regime. The others are assuming a more realistic density profile for the dark 
matter component, and avoiding treating the ionized matter as a fluid by going directly to 
a kinetic description, introducing explicitly the dynamical balance between the recombined 
and ionized fractions of ordinary matter. However, as it stands it already shows that there 
are rich magnetohydrodynamical processes occurring as early as Z ~ 100, which must be 
properly understood before our picture of the physics of cosmic magnetic fields is complete. 
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7. Appendix 

The Christoffel symbols needed to calculate the covariant derivatives 

r 8o = i + ^ n k = ^,k 

rio = % (! - i>) r? fc = [J (i - - V>] 7i* 

In the context of the paper, these are further simplified, since ip — 0. 
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